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Abstract 

We study approximations of reflected Ito diffusions on convex subsets D of M. d 
by solutions of stochastic differential equations with penalization terms. We as- 
sume that the diffusion coefficients are merely measurable (possibly discontinuous) 
functions. In the case of Lipschitz continuous coefficients we give the rate of LP 
approximation for every p > 1. We prove that if D is a convex polyhedron then the 
rate is ©((i^f) 1 / 2 ), and in the general case the rate is ©((i^f) 1 / 4 ). 
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1 Introduction 

In the paper we study weak and strong approximations of solutions of d-dimensional 
stochastic differential equations (SDEs) 

X t = x + [ a(s,X s )dW s + [ b(s,X s )ds + K t , t G R + (1.1) 
Jo Jo 

with reflecting boundary condition on a convex domain D. Here xq G D = D U dD, 
X is a reflecting process on D, K is a bounded variation process with variation \K\ 
increasing only, when Xt G dD, W is a (i-dimensional standard Wiener process and 
a : R + x R d -)■ M. d (g> R d , b : R + x R d -> M. d are measurable (possibly discontinuous) 
functions. Suppose that for n G N we are given measurable coefficients a n : M + xR d -> 
M d (g)]R d , b n : M + x R d — >■ R d and a standard Wiener process W n , and assume that there 
exists a solution X n of the following SDE with penalization term 

X" = x + /* tr n (8,X?)dW?+ f ' b n {s,X^)ds-n f \x?-IL(X?))da, t G R + , (1.2) 
Jo Jo Jo 

where H(x) is the projection of x on D. The problem is to find conditions on {ct^}, {b n } 
ensuring convergence of {^ n } to the reflected diffusion X, and secondly, to give the rate 
of such convergence. 
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Reflected diffusions have many applications, for instance in queueing systems, seis- 
mic reliability analysis and finance (see e.g. Asmussen [3], Dupuis and Ramanan [UJ, 
Kree and Soize [ID], Pettersson [18], Shepp and Shiryaev [21]). Therefore, the problem 
of practical approximations of solutions of (jl.ip is very important. Discrete penalization 
schemes based on the approximation of X by solutions of equations with penalization 
term are well known (see e.g. Pettersson [19] . Kanagawa and Saisho [8], Liu [13] . 
Slomihski [26] ). 

Approximation of reflected diffusions via penalization methods was earlier considered 
by Menaldi [16], Menaldi and Robin [17], Lions and Sznitman [IT], Lions, Menaldi and 
Sznitman [12], Storm [22], Saisho and Tanaka [22] and many others. Unfortunately, 
these authors have restricted themselves to the case of Lipschitz continuous coefficients. 

In the present paper we consider measurable coefficients cr n , b n such that 

\\a n (t,x)f + \b n (t,x)\ 2 <C(l + \x\ 2 ), (t, x) € R + x R d , n G N (1.3) 

for some C > 0. To prove convergence of {^ n } to X we first show that under (jl.3p 
the sequence {^ n } is very close to the sequence {II(X n )} and we observe that H(X n ) 
is a solution of some Skorokhod problem (for the definition of the Skorokhod problem 
see Section 2). Next, using a well developed theory of convergence of solutions of the 
Skorokhod problem (see e.g. [3j[20j[25l[26l[29]) we prove our main approximation results. 
Moreover, we are able to strengthen the rate of the convergence of the penalization 
method in the classical case of Lipschitz continuous coefficients a, b. 
The paper is organized as follows. 

In Section 2 we estimate the L p distance between X n and D. Using some new 
estimates of L p -modulus of continuity of Ito's processes from Fischer and Nappo [7] we 
prove that under (jl.3p for every p > 1, T > 0, 

||su P dist(Xf,L>)li P = e>((— )V2), 
t<T n 

where || • [L = (E(-) p ) 1 / p denotes the usual L p norm. We also show that {AT n } is tight 
in C(M + ,lR rf ) and its weak limit point solve the Skorokhod problem. 

Section 3 contains our main results concerning weak and strong approximations of 
solutions of (jl.ip . We consider the set of conditions on coefficients from the paper by 
Rozkosz and Slomihski [2DJ on stability of solutions of stochastic differential equations 
with reflecting boundary. Roughly speaking, we assume that {a n }, {b n } satisfy (jl.3p and 
{(det ffn^n) -1 } is locally uniformly integrable on some set (<r* denotes the matrix adjoint 
to a n ). Then we show that if {a n }, {b n } tend to a, b a.e. on the set mentioned above 
and uniformly on its completion, then X n — >x> X, where X denotes a unique weak 
solution of (jl.ip . Under the additional assumptions that W n — >-p W and that (jl.ip is 
pathwise unique we show that X n — >-p X. Thus, we generalize earlier approximation 
results to equations with possibly discontinuous and nonelliptic diffusion coefficients 
and discontinuous drift coefficients. 

Section 4 is devoted to the classical case, where all coefficients are fixed Lipschitz 
continuous functions with respect to x and all stochastic integrals are driven by the 
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same Wiener process, i.e. a n = a, b n = b, and W n = W, n G N, and there is L > such 
that 

\\a(t,x) - (7{t,y)\\ 2 + \b(t,x) - b(t,y)\ 2 < L\x - y\ 2 , (t,x)£R + xR d . (1.4) 
In this case, if D is a convex polyhedron, we prove that for every p > 1, T > 0, 

iisup|xr-^iii P = o((— ) 1/2 )- 

For arbitrary convex domain we prove that for every p > 1, T > 0, 

\\ S up\X?-X t \\\ p = 0(( 1 ™)^). 
t<T n 

Thus, we strengthen earlier results on the subject proved by Menaldi [16J. 

In the sequel we use the following notation. R + = [0, oo), C(IR + ,]R rf ) is the space of 
continuous functions x : M. + — > M. d equipped with the topology of uniform convergence 
on compact subsets of M + . For every x G C(M + ,M rf ), S > 0, T > we set ojg(x,T) = 
sup{|x t - x s \; s,t G [Q,T], \s - t\ < 5}. R d ® R d is the set of (d x d)- matrices. The 
abbreviation a.e. means "almost everywhere" with respect to the Lebesgue measure, 
" — >T> n , " — >v" denote convergence in law and in probability, respectively. 



2 General results 

Let Dbea nonempty convex domain in R d and let M x denote the set of inward normal 
unit vectors at x G 3D. Note that n G Af x if and only if < y — x,n >> for every 
y e D (see e.g. [13 [22]). Moreover, if dist(x,Z)) > 0, then 

\Ii{x)-x\ GA/n ^- 

Let Y be an {J-t} - adapted process with continuous trajectories. We will say that a pair 
(X, K) of {Tt} - adapted processes is a solution of the Skorokhod problem associated 
with Y if 

X = Y + K, 
X is D- valued, 

K is a process with locally bounded variation such that Kq = and 

K t = I n s d\K\ s , \K\ t = f l {Xs EdD}d\K\ s , t G R + , 
Jo Jo 
where n s G Mx s if X s G dD. 

It is well known that for every process Y with continuous trajectories there exists a 
unique solution (X, K) of the Skorokhod problem associated with Y (see e.g. [1] or |15j . 
where a more general case of cadlag processes is considered) . The theory of convergence 
of solutions of the Skorokhod problem is well developed (see e.g. [31 [221 ESI ESI 122] )■ 
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Unfortunately, solutions of (jl,2p are not solutions of the Skorokhod problem and the 
problem of their convergence is more delicate. 

Suppose that we are given a filtered probability space (Q n , J 711 , {'FJ 1 }, P n ) satisfying 
the usual conditions and a d - dimensional {J 7 "} -Wiener process W n , n G N. Let {X n } 
denote the sequence of solutions of (jl.2p . In the present paper we will use the simple 
fact that under (11, 3h there exists a sequence of solutions of the Skorokhod problem very 
close to the sequence {X™}. Observe that we can rewrite (|1.2p into the form 

n(X t n ) = Y t n - n [\x™ - U(X^))ds, t G R+ 
Jo 

where Y t n = x - X? + IL(X?) + Jj a n (s, X™) dW^ + jj b n (s, X") ds, t G R+, n G N. 
Since H(X n ) G £>, |i^ ra | increases only when U(X n ) t G S.D and 

K? = n jT |n™I^!| |nTO " I*" = ^ n s i G R+, 

it is clear that (II(X n ), K n ) is a solution of the Skorokhod problem associated with Y n , 
n G N. One can also observe that 

\X? - n(Xf)| = dist(X™, D), t G M + , n G N. 



Theorem 2.1 Assume that $1.3\) is satisfied. 
(i) For every p > 1, T > i/iere is C > suc/i £/iai 

||supdist(X t n ,5)l|L < C(— ) 1/2 , n G N. 
t<T « 

(ii) {(X n , K n )} ne N «s in C(R + ,IR M ) and its every weak limit point (X,K) is a 
solution of the Skorokhod problem. 

Proof, (i) Fix T > 0. First observe that by [14, Corollary 2.4] and Gronwall's lemma, 

sup£sup|X 4 n | p < +oo (2.1) 

n t<T 

for every p > 1. By the above and estimates from Fischer and Nappo Theorem 1] 
(see also [18, Lemma 4.4] and [26, Lemma A4]) for every p > 1 there is C > such that 

\\u 1/n (Y n ,T)\\ p <C(— ) 1 / 2 , nGN, (2.2) 

' n 

where Y t n = x + J * ct„(s, X") dW? + jj b n {s, X s n ) ds, t G M + , nGN. 

Fix n G iV and A; = 0, 1, [nT] — 1. Clearly, X n is a solution of the equation 

X n k , n+S = X n kln + Y k n /n+s - Y k n /n - n J\xi /n+u - U(X- /n+ Jdu, s G [0, 1/n] (2.3) 
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on the interval [k/n, (k + l)/n]. It is also clear that there exists a unique solution of the 
equation 

X? = K/n ~ n f{K ~ K{K))du, s G [0, 1/n]. (2.4) 
Jo 

One can easily check that X% = U(X^ /n ) + - U(X^ /n ) e - ns , s G [0,1/n], which 

implies that 

l*i% - n W/n)l < \ty n ~ n(^/JI = \K,n - n^/Jle" 1 . (2.5) 

Subtracting (12. 4p from (12 . 3|) we see that 

X n k/n+s -X: = Y k n /n+ -Y k n /n -n J\x n k/n+u -X:-U(X n k/n+u +U(X n H ))du, s G [0, 1/n], 
hence that 

\X n k/n+s - X?\ < u 1/n (Y n ,T) + 2nj' \X n k/n+u - X^\du, s G [0, 1/n], 

because II : M rf — > D is Lipschitz continuous with the constant equal to 1. Applying 
Gronwall's lemma we conclude from the above that 

\^ k+l)/n -Xl /n \<e 2 u l/n (Y\T). (2.6) 

Setting A% = |X» - n(X™,| and using ([23]), ((2SD we have 

ai +1 < \x^ k+l)/n - n{xi /n \ < \x^ k+l)/n - xfy n \ + \xfy n - n(X™ /n )\ 

< e 2 u 1/n (Y n ,T) + e~ l A n k . 
Since Aq = and A" < e 2 uJi/ n (Y n , T), by induction on k we obtain 



0< max T] " n(X" /n )| < ^^-^(F", T). (2.7) 

Furthermore, for A; = 0, 1, [nT] and s G [0, 1/n] such that fe/n + s < T, 

yn yn \rn -y-n 

^k/n+s ^k/n ~ 1 k/n+s 1 k/n 



-nj\(x n k/n+u - xs /n ) - (U(X^ /n+u ) - U(X- /n ))d 
-ns{X% /n -IL{Xl /n ). 
Hence, by Gronwall's lemma, 

sup |X ( " fe/n+s)AT - X n k/n \ < e 2 (u 1/n (Y n ,T) + \X» /n - U(X^ /n )\), 

se[0,l/n] 



which when combined with (I2.7P gives 



max sup \X? k/n+s)AT -XJt /n \<Cu; 1/n (Y\T) (2.8) 

0<fc<[nT] s€ [ 0]1 / n ] 



for some C > 0. Of course (|2T7j) . (l2T8j) and (JO} imply (i). 

(ii) By (jl.3p . (|'2.1|) and the well known Aldous criterion (see e.g. [2]), 

{y n } is tight in CQR + ,lR d ). 

Moreover, {? n } satisfies the so called UT condition (see e.g. [251 126| ) and hence its 
every weak limit point is a semimartingale. Due to part (i), the sequence {Y n } is also 
tight in C(M + ,M d ). Assume that Y (n "> — y v Y in C(K+,M d ) along some subsequence. 
By [221 Corollary A3], 

(X^ n \K^)^(X,K) mC(R + ,R 2d ), 

where (X,K) is a unique solution of the Skorkhod problem associated with a semi- 
martingale Y. □ 

Remark 2.2 Under (Ol) . 

sup^liT 1 ^ < +oo (2.9) 

n 

for every p > 1, T > 0. This follows from (|2.ip and Theorem 2.5]. 
3 Approximations of weak and strong solutions 

We say that the SDE (jl.ip has a strong solution if there exists a pair (X, K) of {J~t}- 
adapted processes satisfying (jl.ip and such that (A", iT) is a solution of the Skorokhod 
problem associated with 

Y t = x + [ a(s,X s )dW s + f b(s,X s )ds, tel + . (3.1) 
Jo Jo 

Recall also that the SDE (jl.ip is said to have a weak solution if there exists a probability 
space (f2, T, P), an {Ft}- adapted Wiener process W and a pair of {Ft}- adapted 
processes (X, K) saisfying (|l.ip with W instead of W. 

The following set of general conditions was introduced in Rozkosz and Slomihski 
[20j . We say that condition (H) is satisfied if for some closed subsets H , H\ of M + x M. d 
such that Hi C H , 

• V e >o {(det <7 n <T*) _1 } nG N is uniformly integrable on each bounded subset of H c (e), 

• <r n — y it, b n — y b a.e. on H c = R + x R d \ H, 

• for every (t,x) G Hi (for every (t,x) € i?), 

cr n (t,a; n ) ->• a(t,x), b n (t, x n ) ->• x) 
for all {x n } such that x n — y x (for all {(t,x n )} C -ff such that x n —y x). 

Here H c (e) = R + x R d \ H(e) and JET(e) = H U ili )£ , where #i, e = if J3"i = and 
i!i )E = {z£ R + x R d : infyg^ \z — y\ < e}, otherwise. 
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Theorem 3.1 Assume that hi. 3^ and (H) are satisfied. 
(i) If the SDE U.l\) has a unique weak solution X then 

X n ^X inC(R + ,M. d ). 

(ii) IfW n — > V W in C(R+,R d ) and the SDE [Ll\) is pathwise unique then 

X n ^X inC(R + ,R d ), 

where X is a unique strong solution of \1.1\) . 

Proof. We use notations from the proof of Theorem 12.11 

(i) Our method of proof will be adaptation of the proof of [201 Theorem 2.2]. Since 
K n is a bounded variation process, one can observe that Krylov's inequality used in [20|, 
Theorem 5.1] is still in force, i.e. there exists a constant C depending only on d, R and 
t such that for every non-negative measurable / : R + xR^ R + , 

E / f{a,X?)ds < C||(deta n O~ 1/(d+1) /ll L d+1 ([o,*]xiJ(o,fl)), (3-2) 
Jo 

where t£ = inf{t : \X t \V\K n \ t > R}, B(0,R) = {y G R d : \y\ < R}. By TheoremO^), 
{(X n ,W n )} is tight in C(R + ,R 2d ) and we may assume that (X^,W^) — >v{X,W) 
in C(R + ,R 2d ) along some subsequence, where W is a Wiener process with respect to 
the natural filtration T x,w . By f|3.2|) and arguments from the proof of [20, Theorem 
2.2], 

(X( n \Y( n \W (n) )^(X,Y,W) in C(R+,R M ), 

where Y t = x + Jj a(s, X s )dW s + jj b(s, X s )ds, t € R+. Since Y n = Y n - X n + U(X n ), 
it follows from Theorem I2.1f ii) that 

(X( n \K( n \Y^)-^(X,Y,K) in C(R + ,R M ), 

where (X, K) is a solution of the Skorokhod problem associated with Y . Hence (X, K) 
is a weak solution of (jl.lB and the result follows due to weak uniqness of (II. ip . 

(ii) By using arguments from Gyongy and Krylov [9], to prove that {X™} converges 
in probability it is sufficient to show that from any subsequences (I) C (n), (m) C 
(n) it is possible to choose further subsequences (4) C (/)>( m fc) C (m) such that 
(X^),x( mfc )) — in C(R + ,R 2d ), where I is a process with continous trajec- 
tories. From Theorem I2.1f ii) we deduce that 

{(X«,x( m ),H/^,^ m ))} is tight in C(R+,R d ). 

Therefore, we can choose subsequences C (I), (rn^) C (m) such that 

( X ^\x {mk \ W^, W {mk) ) —>(*', W 7 , m C(R+, R 4d ), 
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where X', X" are processes with continuous trajectories and TV is a Wiener process with 
respect to the natural filtration T x ,x • . In view of part (i), the processes X',X" are 
solutions of (jl.ip with W in place of W. Since (jl.ip is pathwise unique, X' = X" , and 
consequently {X n } converges in probability in C(R+,M 4d ) to some continuous process 
X. Hence (X n , W n ) — >v(X, W), so using once again the pathwise uniqueness property 
of (II. ip shows that X is a unique strong solution of (jl.lj) . □ 

Remark 3.2 Prom |20| it follows that in fact in part (i) of the above theorem the 
assumption that a n —¥ a a.e. on H c may be replaced by a weaker assumption that 
ff„(j* — > aa* a.e. on H c . 

Remark 3.3 There are important examples of equations of the form (jl.ip with discon- 
tinuous coefficients having unique weak or strong solutions. For instance, in Schmidt 
|21j it is shown that if d = 1, D = (t*i, rg), 6 = and a is purely function of x, then (jl.lj) 
has a unique weak solution for every starting point xq € D if and only if the set M of all 
x € -D such that /^ n( y a~ 2 (y) dy = +oo for every open neighborhood C/^ of x is equal 
to the set N of zeros of a. In multidimensional case it is known, that a solution of (jl.ip 
is unique in law if (jl .3p is satisfied with c n ^ b n replaced by a, b, the coefficient aa* is 
bounded, continuous and uniformly elliptic, and dD is regular (see Stroock and Varad- 
han [28] for more details). Recently Semrau |23j considered the classical case d = 1, 
D = R + with coefficients a, b depending only on x. She has shown that if a, b, satisfy 
(|1.3p . a is uniformly positive and (<r(x) — <r(y)) 2 < \f(x) — f(y)\, x,y £ M + for some 
bounded increasing function / : M + — > R then there exists a unique strong solution of 
(jl.ip . Some weaker results on pathwise uniqueness can be found in the earlier paper by 
Zhang [30]. 

Since in condition (H) we do not require continuity of the limit coefficients a and b, 
Theorem 13.11 is a useful tool for practical approximations of solutions of the equations 
mentioned above. 

4 Rate of convergence in the case of Lipschitz continuous coefficients 

In this section we assume that a n = a, b n = b, n S N, where a, b are Lipschitz continuous 
functions with respect to x, i.e. satisfy (jl.4p . We also assume that all SDEs with 
penalization term are driven by a fixed {J-j}- Wiener process W. In particular, this 
means that X n is a solution of the equation 

X t n = x + I a{s, X s n ) dW s + f b(s, X^) ds-n f * (X? - U(X?))ds, t G M+. (4.1) 
Jo Jo Jo 

Tanaka [29] has shown that in the case of Lipschitz continous coefficients there exists a 
unique strong solution (X, K) of (jl.ip . Moreover, from [26, Theorem 2.2] and Gronwall's 
lemma it follows that 

EswY>\X t \ p < +oo and E\K\ P T < +oo (4.2) 

t<T 

for every p > 1, T > 0. 
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Theorem 4.1 Assume that \1.3\) and jl-4\ ) are satisfied. Let X n satisfy \4- 1\) > n£N. 
For every p € N, T > there is C > such that 

(i) if D is a convex polyhedron then 

||sup|X™ -Xt\\\ p < C(— ) 1/2 , n G N, 

(ii) «/ D is a general convex domain then 

\\ S vv\X?-X t \\\ p <C(— ) 1/4 , n G N, 

where X is a unique strong solution of M. 1}) 

Proof. Fix T > 0. Without loss of generality we may assume that p > 2. 

(i) By Theorem 2.2 from Dupuis and Ishi [5] there exists C > such that 

sup |n(X") - Z s | < Csup \Y s n - Y s \ 

s<t s<t 

< C(sup |n(X s n ) - X™\ + sup \Y s n - Y s \) (4.3) 

s<t s<t 

for every t < T, where Y s n = x + / Q s <t(«, XJ)dW u + J S 6(«, Y s n = Y s n + 

n(X^) - X™, s < T, n e N. Therefore, by Theorem O^i), Burkholder-Davis-Gundy 
and Schwarz's inequalities, 

£sup \X^ - X s \ p < Const((— f' 2 + £sup \Y s n - Y s \ p 

S<t n 3<t 

< Const {{—Y> 2 + E / | \o-(s,X™) -a(s,X s Wds 
n Jo 

+E [ \b(s,X^-b(s,X s )\Pds) 
Jo 

for every t <T. By the above and (|1.4p . 

Esup \X? - X s \p < Const ((— )P/ 2 + f Esup \X% - X u \ p ds) 
s<t n Jo u<s 

for every t < T, so (i) follows by Gronwall's lemma. 

(ii) If D is a general convex domain then by Lemma 2.2 in Tanaka [29], 

|n(xp) - x 4 | 2 < |F™ - y| 2 + 2 /V/ 1 - >* - + y,) d(K? - k„)\ 

Jo 

< ConstflnpC?) - X t n \ 2 + \Y t n - Y t \ 2 + sup |II(X t n ) - X t n \{K n \ T + \K\ T ) 

t<T 

+ | A?" - Y t - Y? + y s ) d(K- - K s )\ (4.4) 
•/ 
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for every t <T. Since by the integration by parts formula, 

f\Y t n -Y t -Y s n + Y s )d(K«-K s ) 
Jo 

= J\x- - X s )d(Y s n ~ Y s ) + \([Y n - Y] t - \Y t n - Y t \ 2 ) 

(here \Y n — Y] denotes the quadratic variation of Y n — Y), it follows from Theorem 
E3P) and gaj) that 

Esup \X? -X s \ p < Const((— ) p / 2 + £sup \Y s n - Y s \ p + E([Y n - Y] t ) p ' 2 

S<t n 8<t 

+E(snp \U(X t n ) - X?\) p/2 (\K n \ t + \K\ t y' 2 

t<T 

+Esup\ f\x^-X s )d(Y s n -Y s )\P/ 2 ). 

s<t JO 

By Schwarz's inequality, Theorem I2.1f i). (|2.9p and (|4.2p . 

£(sup |npf") - X?\)P/ 2 (\K n \ T + \K\ T y' 2 < Const((— ) p / 4 ). 

t<T n 

Since Y n — Y is a continuous semimartingale with a martingale part M n = J Q a(s, X™) — 
cr(s, X s )dW s and a bounded variation part V n = j Q b(s,Xg) — b(s,X s )ds, using Burk- 
holder-Davis-Gundy and Schwarz's inequalities we get 

£sup | [\x? - X s )d(Y s n - Y s )\ p ' 2 < Const(£( f \X? - X s \ 2 d[M n ] s ) p ' A 

s<t Jo Jo 

+E(su P \xz - x s \\v n \ t y' 2 ) 



s<t 

< Const(£sup |X S " - X s \P)V 2 (E([M n ] t y/ 2 + {\V n \ t f) l l 2 . 

s<t 

Observing that [Y n — Y] = [M n ] and using the elementary inequality 2ab < e 2 a 2 + (b/e) 2 
with some sufficiently small e we deduce from the above that 

Esup |X s n - X s \ p < Const((— ) p / 4 + Esup \Y s n - Y s \ p 

s<t n s<t 

+E([M n ] t ) p / 2 + E(\V n \ t ) p )) 

< Const ((— ) p/4 + £ / \\a(s,X^)-a(s,X s Wds 

n Jo 

+E [ \b(s,X^)-b(s,X s )\ p ds) 
Jo 

< Const ((—f/ 4 + / Esup\X% ~X u \ p ds) 

n Jo u<s 

for every t <T. Using Gronwall's lemma completes the proof. □ 
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Remark 4.2 In the case of bounded convex domains and bounded Lipschitz continuous 
coefficients a, b the problem of L p approximation of solutions of (|l.ip by sequences of 
solutions of (|4,ip was considered earlier in Menaldi [16]. In particular, in (161 Theorem 
3.1] it is proved that for every p > 1 and T > 0, || sup t<T [AT™ — At||| p — >• 0. From the 
proof of \16\ Theorem 3.1] one can also deduce that 

V 5>0 ||sup|Xr-X t j|| p = 0((i) 1/4 ^)- (4-5) 

t<T n 

In fact, in [16^ Remark 3.1] a better rate is stated. However, R. Pettersson has observed 
that there is a gap in the proof of Theorem 3.1] (in the first line on page 741 p should 
be replaced by 2p). Using Menaldi's calculations and taking into account Pettersson's 
remark one can only prove (|4.5p . It is also worth pointing out that Menaldi's method of 
proof of (|4.5p is completely different from our method based on estimates of L p -modulus 
of continuity for Ito processes. 
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